iBiw 




.ota 


Copy 




a ***i**~ 

[W -' 


>U2uC* 

By authwsi 


i>iU 50 UNCLASSIFIED^ 



RESEARCH MEMORANDUM 


PRELIMINARY EMPIRICAL DESIGN REQUIREMENTS FOR TEE 
PREVENTION OF TUMBLING OF AIRPLANES 


HAVING NO HORIZONTAL 


taEPB Eeferenc 


By Robert L. Bryant • JJ J 1 111 ~ 

POI TO BE X&KET PBCK THH EOQ#t 

Langley Aeronautical Laboratory 
Langley Air Force Base, Va. 

SSIFICATION CANCELLED^ CLASSmcATO eHAHOSD 


to 


ttda document contains classified 

me a nin g of the Eaptaage Act, U3C SCtM Lt» ref l ation of Its consents (fa Jy ' * l t 

manner to an unauthorised per son Is proBfbCted by law, w 

Information so cl a ssifi ed nay be Imparted only to persona la the military ud naval services of the Craltad 
States, appropriate civilian officers tod employees of the Federal O mnmM flt who have a legitimate Into root 
therein, and to United States citizen* of kncwn loyalty and discretion wha of necessity must he tnJbrmed thereof. 

NATIONAL ADVISORY COMMITTEE 

FOR AERONAUTICS UNCUSS/Flf 



WASHINGTON ^ w A UJ4 r '-KY 
October 11, 195CM*"^^ ! ^ 





NACA EM L50H23 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

PRELIMINARY EMPIRICAL DESIGN REQUIREMENTS FOR THE 
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An investigation has been made of the design characteristics and 
loadings that are conducive to the tumbling of airplanes that have no 
horizontal tails. Preliminary empirical design requirements based on 
model tests of 18 different configurations are presented. A brief 
explanation of the phenomenon of tumbling is appended. 


INTRODUCTION 


The summary of results of tumbling investigations of several dynamic 
models to determine their tumbling characteristics (reference l) indi- 
cated that- tumbling may occur only for airplanes having no horizontal 
tails and that the degree of. static longitudinal stability was indica- 
tive of the proneness of an airplane to tumble . From analysis of the 
data of reference 1 and of additional data obtained in the present 
investigation certain empirical design requirements for the prevention 
of tumbling of horizontal tailless airplanes were indicated and are 
presented herein. Results were also obtained in the present investi- 
gation from tumbling tests of five flat plates. Tumble theory advanced 
by Dupliech in reference 2 is also discussed and appended herein. 

Although further research will be necessary for complete solution of 
the problem, it appeared desirable to present an empirical criterion at 
this time based on existing knowledge . 


SYMBOLS 


All angles and vectorial quantities are shown in figure 1. 

I distance measured along the horizontal earth axis, feet 
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distance measured along the vertical earth axis, feet 
angle of attach, degrees 

flight-path angle relative to g earth axis, degrees 

angular displacement in pitch with respect to the earth g 
axis, degrees 

average rate of pitch during any cycle of the tumble after 
equilibrium has been obtained, radians per second 

instantaneous rate of pitch, radians per second 

instantaneous resultant velocity, feet per second 

p 

acceleration tangent to flight path, feet per second 
acceleration normal, to flight path, feet per second^ 
time, seconds 
mass, slugs 
weight, pounds 

moment of inertia about the lateral axis, slug-feet^ 

p 

wing area, feet 
wing span, feet 

wing chord at any station, feet 
mean aerodynamic chord, feet 

distance from center of gravity to centroid of plan-form area, 
feet 

air density, slug per cubic foot 

relative-density coefficient based on span 

resultant aerodynamic force, pounds 
resultant inertia force, pounds 
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L 

D 

M 

c L 


lift, pounds 
drag, pounds 

pitching moment, foot-pounds 


lift coefficient 





J N 


pitching-moment coefficient 


Pitching moment'' 


ipV 2 sc 


normal-force coefficient 


( ! 

( formal force\ 

lpv% ) 


longitudinal-force coefficient 


f Longitudinal force 


lpV 2 S 


: ) 


Parenthetic subscripts: 

( ) a due to angle of attack 

( )* c due to time rate of change of angle of attack 

( )g due to pitching velocity 


APPARATUS AMD TESTS 
Models 


The models vere scaled from full-scale airplanes with the exception 
of two research models G and I and five flat plates, models M, 0, F, Q, 
and R, which were built to an arbitrary scale. The dimensional charac- 
teristics, mass characteristics, pitching-moment characteristics, and 
plan views of the models are given in tables I and II in terms of full- 
scale values. The models were made of balsa and ballasted with lead 
weights to obtain dynamic similarity to their actual or assumed full- 
scale counterparts at some arbitrary altitude as indicated on table II. 
The assumed full-scale counterparts of the research models and flat 
plates were such that their spans were within the range of the other 
models tested. 
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Wind Tunnel and Testing Technique 

The. tumble tests were conducted in the Langley 20-foot free-spinning 
tunnel, a vertical wind t unne l of dodec agonal cross section capable of 
airspeeds up to approximately 60 miles per hour and is similar to the 
former Langley 15-foot free-spinning tunnel as described in reference 3* 
For the tests, two methods of launching the models were employed: the 

model was released from a nose-up attitude to simulate a whip stall in 
order to determine whether the model would start tumbling of its own 
accord, or the model was launched with initial pitching rotation in order 
to determine whether the model would continue to tumble once the tumbling 
motion had been started. The simulated whip stall was obtained by holding 
the model in the air stream with its nose up and simply letting go of the 
model. When launched with initial pitching rotation, the model was held 
in the air stream and forced to rotate by applying a torque about the 
lateral axis as it was launched. 

The models were generally tested with three elevator deflections, 
full up, neutral, and full down with both positive and negative initial 
pitching rotation. Whip-stall tests were generally made only for those 
conditions for which the model tumbled when given forced rotation. It 
was considered that the model would tumble for a particular center-of- 
gravity location if a tumble was obtained by any method of launching for 
any elevator position. 


Parameters Considered for Empirical Criterion 

The results of the investigations summarized in reference 1 indi- 
cated that aspect ratio, the square of the radius of gyration in percent 
of the span squared and the degree of static longitudinal stability 
appeared to be important parameters with regard to tumbling. When in 
the analysis of available data and data obtained in the present investi- 
gation no way of using these parameters was found that would fit into the 
scheme of separating the tumbles from no tumbles, it was decided to look 
for parameters whose effect on the tumble were most obvious. Variation 
of the center-of -gravity location has a decided effect on tumble results. 
The center-of-gravity location determines the degree of static longi- 
tudinal stability which, in conjunction with the moment of inertia, 
deteimines the initial rate of pitch. Since the rate of pitch, which, 
in the case of tumbling is equal to the rate of change of angle of attack, 
determines the amount of hysteresis in lift (see appendix), it was 
decided to consider the center-of-gravity location in percent of the 
mean aerodynamic chord as a basic parameter. The effect of the damping 
moment is quite obvious inasmuch as conventional-type models (with 
horizontal tails) which have a large degree of damping will not tumble 
even when the center-of-gravity location is such that extreme static 
instability exists. If the damping moment always exceeds the moment 
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due to hysteresis in lift for any rate of change of a ng le of attack, the 
t umb ling motion can never get started. The usual expression for wing 
damping which considers only relatively small rates of pitch may not he 
applicable for tumbling motions in which the rates of pitch are appreci- 
ably larger. It was felt that some measure of the damping would be indi- 
cated by the square of the distance between the center of gravity and 
the centroid of the plan-form area. Damping, in addition to being a 
function of the rate of pitch, is also a function of the rate of descent; 
the rates of pitch and descent are respectively dependent upon moment of 
inertia and mass, and therefore it was felt that the ratio of mass to 
moment of inertia, or radius of gyration sauared should also appear in 

the damping parameter. The expression y y — or y— — was thus chosen 
to represent damping. 


RESULTS 


The results of the tumbling investigation for the 11 specific models 
and the 2 research models are presented in figures 2 and 3* Figure 2 
presents the results of tests of models having aspect ratios greater 
than 3 and varying up to 11, whereas figure 3 presents the results of 
tests of models having aspect ratios of 3 or less. The results of the 
tumbling tests of the flat plates are shown in figure 4. The coordinates 
on the figures are the two parameters previously discussed, center-of- 
gravity position and mh^/ly. Results for both forced rotations and for 
whip stalls are shown. 

Interpretation of model tumbling results as applied to corresponding 
full-scale airplanes regarding the possibility of tumbling is as follows: 
If a model does not tumble when launched with forced rotation, the corre- 
sponding airplane will not tumble; if the model tumbles when launched 
with forced rotation the corresponding airplane may tumble but there will 
be greater likelihood of tumbling if the model tumbles when launched from 
a whip stall. 


DISCUSSION 


Where data permitted, two separation lines have been dr?wn in fig- 
ures 2 to 4 to separate the data for the models which tumbled when given 
forced rotation from the data for the models which tumbled from a whip 
stall and from the data for the models which did not tumble. In addi- 
tion, a line has been drawn, where possible, below which tumbles from a 
whip stall appear unlikely. The dashed lines are conjectural and have 
been drawn on the assumption that the separation lines are a family of 
curves as some function of aspect ratio. 
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Tumbling tendencies of nine models whose aspect ratios were greater 
than 3 and ranging up to 11 are indicated in figure 2. At a given value 
of the parameter mh^/ly, the models would not tumble even when launched 
with forced rotation for relatively forward center- of -gravity positions. 

As the center of gravity was moved rearward, tumbles were obtained first 
only from forced rotation, then from forced rotation and occasionally 
from whip stalls, and finally it appeared likely that tumbles would 
always be obtained even from whip stalls for rearward center-of -gravity 
positions. At a forward center-of-gravity position, increasing the 
value of mh^/iy tended to decrease the tumbling tendency from forced 
rotation, whereas for a rearward center-of-gravity position the tendency 
to tumble from whip stalls was increased for an increase in mh2/ly. 

This difference in effect of mh^/ly may be partially attributable to 
a difference in rotational inertia. Large radius of gyration (high 
inertia with low mass) would tend to prevent the tumbling motion from a 
whip stall, whereas once the tumbling motion has been obtained from 
forced rotation, it would aid in sustaining the motion. 

These results indicate that the center of gravity need not be 
nearly so far forward to avoid tumbles from whip stalls, as it must be to 
avoid tumbles from forced rotation at low and moderate values of mh2/iy. 
Thus, although no comparison between model and airplane tumbling results 
is available, it appears that unless an airplane encounters an extrinsic 
force in flight (similar to that given a model when launched with forced 
rotation, which is believed to be a rather remote possibility) the center 
of gravity indicated on the figure as necessary to prevent tumbling may 
be conservative. Further, references 4 and 5 indicate that the increase 
in lift due to pitching velocity decreases with either an increase in 
Reynolds number or Mach number. There- is one instance (data unpublished) 
where a body of near conical shape rotated in a free stream about an 
axis normal to the axis of symmetry but did not rotate at 1 when either 
the Reynolds number or Mach number was increased. In spite of the fact 
that a cone is not a very efficient lifting device, these results indi- 
cate that some scale effect may be expected and that the full-scale air- 
plane may be less prone to tumbling than the model. 

Figure 3 presents the results of tumbling investigations for four 
models whose aspect ratios were 3 or less. For these models, tests were 
made only from forced rotation, and again the models would not tumble 
for relatively forward center-of-gravity positions but would tumble for 
rearward positions. 

Comparisons of results in figures 2 and 3 indicates that, for the 
models of low aspect ratios, further rearward center-of-gravity positions 
were peimissible than for models with large aspect ratio before tumbling 
tendencies were obtained. These results appear to be consistent with those 
in reference 2. The reason for the difference due to aspect ratio was not 
apparent • . - ......... — ‘ - 
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Results of tumbling tests of five flat plates axe presented in fig- 
ure 4. These results show separation between conditions of no tumble 
and those of tumbles similar to those obtained for the models. The 
center-of-gravity positions for which no tumbles were obtained were 
farther rearward, however, than were those for the specific models of, 
similar aspect ratios (fig. 2). As previously noted, however, a specific 
model was said to be capable of tumbling if tumbles occurred for anfy 
elevator position which includes elevator set to aid the tumbling Rota- 
tion. Because the tendency to tumble is greater when the elevators are 
deflected in the direction of the forced rotation than when the eleva- 
tors are neutral, and inasmuch as the flat plates had no elevator, (the 
results of tests of the flat plates (with essentially neutral elevator) 
may therefore be somewhat optimistic with regards to the center-of- 
gravity position required to prevent tumbling. The differences between 
the results of tests of flat plates and specific models may also be in 
part due to the baBic physical differences between flat plates and 
actual airfoil sections. As noted in the appendix, the hysteresis in 
lift is considered to be the primary motive force in tumbling motions. 

The results of reference 6 indicate that thin airfoils have less hyster- 
esis in lift than thick airfoil sections. It would appear therefore 
that flat plates and wings of thin airfoils may generally have less 
tendency to tumble than wings of thick airfoil sections. 

The results obtained with the flat plates indicate that in the 
region in which tumbling occurred only by forced rotation, quite a 
number of no-tumbling conditions were also obtained. This is probably 
due to the fact that the parameters used for the current empirical 
criterion to obtain separation between tumbling and no-tumbling condi- 
tions do not completely dictate this separation and that other factors 
may have some influence. 

Some of the models tested tended to assume motions other than 
tumbling motion when launched for tumble tests. Such motions as rolling 
or cart wheeling (rotation about the Z-axis) were experienced. The 
results of all tumbling tests indicate that the models tend to tumble 
or rotate about the axis of least moment of inertia and/or aerodynamic 
damping. If, for example, the moments of inertia about the X and Y body 
axes were the same, some models may tend to tumble and others to roll, 
depending on the relative damping. Nonuniform gyrations about all axes 
may also be possible. No investigations have been made of these specific 
or different motions. As an example, however, model H, has been observed 
to roll about the X body axis rather than to tumble. This rolling motion 
could be stopped by aileron movement against the rolling motion. It Is 
probable that other nonuniform motions which have been observed could be 
controlled by proper control manipulation Inasmuch as the models generally 
assumed conditions of unstalled flight during the course of the motion. 
Such nonuniform motions generally occurred when the moment of inertia 
about the Y-axis was larger than that about the X-axis. 
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In summary, the separation curves for the specific models (figs. 2 
and 3 ) have been redrawn in figure 5 showing an empirical criterion for 
the prevention of tumbling. Design of airplanes without horizontal tails 
in accordance with this criterion should reasonably assure that the air- 
plane will not tumble provided it is within the range of geometric and, 
mass characteristics of the designs presented herein. It appears desir- 
able to make the design fall- below the lowest criterion curve (fig. 5) j 
but, should this be impracticable, the design should fall below the 
separation curve below which tumbles from whip stalls are not likely. 
Center-of-gravity positions that fall below this line are generally 
necessary for normal stability. 


CONCLUDING REMARKS 


The results of tumbling investigations of several models in the 
Langley 20 -foot spinning t unn el have been studied with regard to the 
design characteristics and loadings that are conducive to tumbling. A 
preliminary empirical criterion for the prevention of tumbling is pre- 
sented which indicates center-of-gravity positions and values of other 
parameters which should prevent tumbling. It appears that further 
research will be required in order to understand completely the problem 
of tumbling. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX 

THEORY OF TUMBLING MOTION 


Duple ich has shown in reference 2 that tumbling is caused by a 
phenomenon associated with a discontinuity in the aerodyn ami c forces 
resulting from a rapid change in angle of attack through, the stall. 

Figure 6 , taken from reference 6 , shows that this phenomenon commonly 
known as hysteresis can. produce large variations in the aerodynamic 
forces due to rate of change of angle of attack. When the angle of 
attack is being varied so that a wing is going from an unstalled flow 
to stalled flow, the flow is altered in such a mann er that greater forces 
result at any given angle of attack than would for a static test at the 
same angle of attack. Conversely, lesser forces result when a wing is 
being unstalled. 

The mechanics of the tumbling motion may be described by considering 
the energy produced by the aerodynamic forces. Shown in figure 7 are 
hypothetical curves of pitching moments against angle of attack which 
are believed to be reasonable and representative of a flat plate in 
tumbling equilibrium with the center of gravity at the 50 percent chord. 
The areas outlined by the moment curves represent the energy resulting 
from the various moments. The flat plate Is used to illustrate the 
mechanics of tumbling inasmuch as it Is bisymmetrical and therefore 
requires only one-half revolution in pitch to complete a cycle of the 
pitching moments. There are three separate aerodynamic pitching moments 
in action during a tumble: the static moment (M) a which initiates 

the pitching motion; the moment due to rate of change of angle of . 
attack (M) ^ which supplies the motive energy that causes the pitching 

motion to accelerate more rapidly; and the moment due to pitch (M)q 
which is the damping moment that causes the pitching motion to reach 
the steady state. At the start of the motion more energy is supplied 
by (M),^ than can be absorbed by (M) q , thus 




(la) 


This Inequality exists until when during the nth cycle of the motion the 
energy supplied by (M)^ is exactly equal to the energy absorbed by (M) q 

and equilibrium Is obtained 
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nnjt+| 

pnir+it 


1 [(M) a + (M)^ + (M) c 

U n n 

J da + / [(M) a + (M) & + (M)J da = 0 

J njt+* 

(lb) 

The kinetic energy of the 
motion or any cycle after 
the total energy supplied 

tumbling motion during the nth cycle of the 
equilibrium has been established is equal to 
in obtaining equilibrium. 


* 

nnir 


I lyq 2 = 

/ [(M) tt + 00 & + (M) q ] da 
u 0 

(lc) 


The energy due to the static moment absorbed by a flat plate or any 
bisymmetrical section during a cycle will be zero if the path of the 
center of gravity is along a straight line. Experiment has shown (refer- 
ences 1 and 2 ) that the deviation of flight path from a straight line is 
small (flight-path angle is nearly constant) and the flight path therefore 
may be considered as a straight line. Equations 1 may then be written 
omitting (M) a . It thus appears that tumbling occurs when the energy 
supplied by the pitching moment due to rate of change of angle of^attack 
exceeds the energy absorbed due to the dancing moment. This condition 
can exist only if the static moment and the moment of inertia are such 
as to permit a rate of change of angle of attack that for the given 
aerodynamic shape will generate sufficient hysteresis in lift. In the 
case of a wing whose static pitching -moment curve is unsymmetrical more 
or less hysteresis in lift must be generated to balance the energy 
resulting from the asymmetry of- the static pitching moment. 

The flight path of a tumbling wing is determined by the balance of 
the instantaneous aerodynamic and inertia forces and moments. (See 
fig. 1.) The equations of motion for the three degrees of freedom are: 


“ 0 - W - j^a + < B >& + < D >q] + [< L >a. + < L >4 + < L )q] 

m - RL) a + (L) 4 + (t) q j sin 7 + [(D) a + (D) 4 + (D)J 


cos 7 


( 2 ) 

(3) 



[(M) a + (M) & + (M)^j 


W 


cos 7 
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Because the variables in equations (2), ( 3 ) > and (4) are nonlinear 
and are interdependent, a general solution of the equations is not 
possible by currently known methods. Because the aerodynamic forces 
and moment are not constant, the velocity of the center of gravity as 
well as the rotation about the center of gravity can not be uniform. 

As previously indicated, however, the resulting flight path will not 

deviate appreciably from a straight line and therefore In 

ax dt 

the case of tumbling, this rotation in pitch ^ is not small, as in 

most stability problems, and exceeds all known values for which data 
have been taken with regard to its effect on lift, drag, and pitching 
moment. It appears that lift, drag, and pitching-moment data must be 
obtained throughout the range of values of pitching velocity peculiar 
to the tumble for any particular airplane in question before a solution 
to equations (2), (3)* and (4) can be obtained. For complete solution 
of the tumble problem more knowledge will be required about the mechanics 
of, and- the parameters which control, hysteresis and damping at values of 
pitching velocity peculiar to the tumble. 
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T ABUS I.- DIMENSIONAL CHARACTERISTICS OP MODELS TESTED 
{Model values are presented in terms of full-scale values] 


Model 

Scale 

Span 

(ft) 

Area ; 
(ft 2 ) 

Aspect 

ratio 

Taper 

ratio 

Sweep l /4 c 
(deg) 

M.A.C. 

(in.) 

L.E. M.A.C. 
aft L.E. root 
(in.) 

D 

l/l 6 

23.30 

427.00. 

1.27 


0.0 

238.00 

10.02 

B 

1/20 

60.00 

490.00 

7.36 

0.248 

21.9 "back 

109.80 

69.70 

C 

1 / 57-33 

172.00 

4020.00 

7.36 

0.248 

21.9 tack 

315.00 

200.00 

D 

1 / 17.8 

54.00 

356.00 

8.20 

0.376 

15.0 fwd 

85.82 

- 29.10 (fwd) 

E 

1 / 17.55 

38.67 

496.00 

3.01 

0.600 

35-0 hack 

157.00 

83.56 

F 

1/20 

51.67 

362.00 

7.37 

0.249 

22.0 back 

94.50 

60.10 

* 

1/20 

44.58 

517-40 

3.84 

1.000 

30.0 hack 

140.00 

77.20 

H 

1/20 

29. 42 

375.00 

2.31 

0.000 

37.6 hack 

203.90 

101.90 

IB 

1/30 

^O.OO 

673.20 

1.34 

1.000 

45.0 back 

272.00 

88.80 

D 

1/16 

39.00 

293.31 

5.19 

0.253 

24.9 hack 

102.30 

49.30 


1/30 

134.00 

1800.00 

10.00 

0.167 

11.2 hack 

190.80 

82.80 

D 

1/60 

290.00 

7920.00 

10.60 

0.200 

0.6 fwd 

375.60 

20.88 

M 

1/20 

26.83 

200.00 

3.60 

0.455 

38-1 back 

93-68 

62.48 

% 

1/20 

50.00 

250.00 

10.00 

1.000 

0.0 

60.00 

0.00 

20 

1/20 

45.00 

225.00 

9.00 

1.000 

0.0 

60.00 

0.00 

2 p 

1/20 

40.00 


8.00 

1.000 

0.0 

60.00 

0.00 

2 Q 

1/20 

35.00 

175.00 

7.00 

1.000 

0.0 

60.00 

0.00 

2 r 

1/20 

40.00 

266.70 

6.00 

1.000 

0.0 

80.00 

0.00 


^Research models. 

p 

Flat plates. 

Others - models of specific airplanes . 
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O fumble (even from whip stall) 
□ Tumble. ( from forced rotation ) 
+ No tumble (even from -forced 
rotation ) 

Letters denote models given 
in table I. 
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Figure 2.- Tumbling tendencies of airplane models having no horizontal 
tail and vhose aspect ratios are between 3 and 11. 
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Figure 3 .- Tumbling tendencies of airplane models having no horizontal 
tail and vhose aspect ratios are between 1 and 3* 
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Figure k.- [fumbling tendencies of flat plates whose aspect ratios are 

"between 6 and 10. 
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Figure 5«” Prel iminar y design requirements for the prevention of 
tumbling of airplanes having no horizontal tail. 
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Figure 6.- Normal- and longitudinal-force coefficients of a Clark YH 
wing whose angle of attack is being varied rapidly at an airspeed 
of 37*5 feet per second. 
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Figure 7«- Hypothetical curves of pitching -moment distribution for a 
flat plate in tumbling equilibrium with the center of gravity at 
the centroid. The areas above and below the zero moment represent 
the energy gained and lost, respectively, during a half -revolution. 





